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Abstract. 



We study the order parameter of a quasi-2D gas of ultracold atoms trapped in an 
optical potential in the presence of controllable disorder. Our results show that disorder 



^ ' drives phase fluctuations without significantly affecting the amplitude of the quasi- 



condensate order parameter. This is evidence that disorder can drive phase fluctuations 
in 2D systems, relevant to the phase-fluctuation mechanism for the superconductor- 
to-insulator phase transition (SIT) in disordered 2D superconductors. 
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1. Introduction 

Although it has been 25 years since the discovery of high-temperature (high-T c ) 
superconductivity in cuprates, the nature of the superconducting phase transition 
surprisingly remains an open question. This is in part due to the effects of disorder 
in two-dimensional (2D) systems: 2D physics is well understood in the context of 
Berezinskii-Kosterlitz-Thouless (BKT) theory [TJ [2], but the addition of disorder can 
make even the phase diagram uncertain [3j. Disorder-driven phase fluctuations in 2D 
systems may produce a Bosonic-type phase transition between a superconducting and 
insulating state in superconducting thin films [U El E] ; there is mounting evidence that 
many phase transitions in high-T c superconductors are indeed Bosonic [7J [8] in nature. 
As it is difficult to access the phase of the order parameter in these systems, other 
mechanisms behind the transition (e.g. pair-breaking) are difficult to rule out. This 
theoretically challenging problem is well suited for quantum simulation. We replace the 
high-T c superconducting system with a well-controlled ultracold Bosonic system, where 
we gain access to the phase of the order parameter through matter- wave interference. 

Two-dimensional systems behave very differently from 3D systems. For instance, 
in 2D, there can be no long- range order. However, this does not preclude the 
existence of a superfluid state for low temperatures, with this superfluid destroyed 
at a critical temperature Tbkt- At non-zero low temperatures, thermal excitations 
exist only in the phase of the order parameter of the superfluid, taking the form of 
long- wavelength phonons and vortices - quantized circular superfluid flow around a 
small spot of non-superfluid. Below Tbkt, vortices are bound into pairs with opposite 
direction of circulation, and the transition to a normal fluid (or gas) occurs after these 
vortex pairs unbind. These phase fluctuations are believed to be the mechanism for 
the BKT transition observed in 2D ultracold atom systems [H [10], as well as the 
temperature-driven SIT in superconducting thin films jU E], and may be relevant 
to the zero-temperature superconducting doping phase transition in high-T c cuprate 
superconductors [7j. In this work, we use laser speckle to controllably add disorder to 
a quasi-2D Boson system of Rb atoms trapped in two planes of an optical lattice. By 
analyzing fringes due to matter-wave interference of the released atoms, we show that 
disorder predominantly drives phase fluctuations, but not fluctuations in the amplitude 
of the quasi-condensate order parameter, expected to be partially superfluid. This work 
complements other studies of disorder in ultracold atom systems [H] in ID [121 IB], 2D 
[H], and 3D p2]. 

2. Experiment 

Our experiment begins with a 3D Bose-Einstein condensate (BEC) of 87 Rb atoms 
confined in a crossed-beam optical dipole trap of frequencies (co x , u y , uj z )/2t: « 
(125, 125,30) Hz. A shallow-angle optical lattice divides this 3D BEC into two nearly 
identical quasi-2D ultracold atom samples confined at two lattice nodes [9] separated 
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Figure 1. Overview of the experimental setup. The atoms are confined in 
a combination of optical trap, shallow-angle optical lattice, and disordering beam. 
Gravity is in the —y direction, while the two planar systems are perpendicular to the 
x-direction. The camera used to take images is situated above the sample. 

by 3.1 /mi, while an additional beam introduces a disordering potential (Fig. [TJ). We 
adjust the temperature of the two samples between 25 and 400 nK by changing the 
final optical trap depth, which simultaneously tunes the total atom number from 15 X 
10 4 to 80 X 10 4 . The lattice provides a confinement frequency in each of the planes of 
u x /2tt = 7.64(4) kHz, which is greater than our largest chemical potentials \ijh ~ 3 
kHz (h Planck's constant) and similar to k^Tjh (ks Boltzmann's constant) at our 
highest investigated temperatures T, indicating that each plane is in the quasi-2D limit. 
The dimensionless 2D interaction constant g ~ 0.21 [16J. Tunneling between planes is 
negligible on the timescale of the experiment, and small occupations in lattice nodes 
outside the central pair are insignificant. 

The disorder in the potential experienced by the atoms is directly proportional to 
the intensity in the disordering beam [TTJ. In order to create spatially small intensity 
features across our ~20x90 /im planes (estimated Thomas-Fermi diameter of the planes 
using the perpendicular lattice frequency [32]), we pass a 777.3 nm laser beam through 
a phase diffuser and focus the beam to a Gaussian waist of 90 /mi on the atoms using 
a high numerical aperture lens (NA = 0.45). The measured Gaussian waists of the 
autocorrelation function of our disorder [18] are 3.2 /tm along the beam axis and 530 
nm transverse to the axis, providing an accurate measure of the (randomly varying) 
average feature size in each direction (see Supplementary Material). Since the disorder- 
inducing beam is perpendicular to the 2D planes, the feature size across the planes is 530 
nm. Additionally, as the spacing between the planes is on the order of the correlation 
length along the beam, we expect that each plane experiences approximately the same 
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disordered potential. 

The 777.3 nm lattice light is turned on with two linear ramps of 500 ms and 100 ms 
to maintain adiabaticity, followed by a 200 ms hold to ensure equilibrium. During the 
first (second) linear ramp, the lattice attains a final potential depth of Uq/H = 6.1±0.3 
(290±20) kHz. The disorder strength is changed linearly from to its final value Af/h 
between and 4.2±0.4 kHz (determined by the final power in the disordering beam) 
during the second ramp. After the hold time, we shut off all of the beams in <200 fis. We 
measure the expanded density distribution of the planes with absorption imaging along 
the ^-direction after 22 ms of expansion and extract information about the distribution 
and interference contrast. 

In the absence of disorder, before shutting off the potentials, the density distribution 
smoothly decreases from a maximum at the center of the cloud to zero at the edge. Using 
a local density approximation (LDA), the local properties of the system are a function 
only of the position-dependent 2D phase space density D(r) = n(r )X 2 dB (T), with n the 
density and Xjb the thermal deBroglie wavelength [HI [19]. At the edges of the cloud, 
where D is low, the gas is thermal. As D increases towards the center of the sample, the 
density fluctuations are suppressed, leading to a pre-superfluid phase which supports 
vortices and can be considered a quasi-condensate - a BEC with a fluctuating phase 
[201 EH I2H [22]. Finally, when D > D c (equivalent to T < T BKT at the local density), 
the system undergoes the BKT phase transition to a superfluid as thermally activated 
vortices pair (as observed in [9]). 

When the potentials are shut off, the clouds expand very quickly along the direction 
of tight confinement (x), producing a nearly interaction-free expansion [231 [201 [23] . I n 
the ^-direction, we see a bimodal distribution (Fig. [2b, d). Atoms in the low-density 
region of the cloud expand ballistically according to the temperature of the sample, 
producing a broad Gaussian background. Atoms in higher density regions near the 
center of the cloud will expand with a characteristic momentum set by the inverse 
position-space coherence, giving a sharp peak in the expansion along z [HI [20]. All 
of our data has -D(O) > D c w 7.4 [251 HB], so there are atoms in the central peak in 
all images. Along x, the two planes completely overlap after expansion and display an 
interference pattern [9l[26] (Fig. [2]). If all of the atoms were 100% Bose-condensed, then 
the contrast of the interference fringes would be 100% [27]. Thermal phase fluctuations 
prevent BEC formation and decrease this visibility through both phonons and vortices, 
even after subtracting off the broad background of non-coherent atoms. The properties 
of the non-disordered trapped quasi-2D Bose gas and its expansion properties have been 
extensively studied E31 EDI EH US [24] . 

We now discuss the effect of disorder on the properties of a quasi-2D Bose gas. 
The presence of disorder can modify both the properties of the trapped gas and the 
properties of the expansion. The density distribution in the trap should take the shape 
of the disorder (more atoms settle into the wells of the disorder), although this effect 
is expected to be screened by repulsive interactions between the atoms [28J. We expect 
that the disorder will also disrupt the phase coherence of the trapped sample. We 
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Figure 2. Sample absorption images. Images after 22 ms expansion show the 
effect of a vortex on the fringe pattern (a), straight fringes (b), and the effect of 
A/h = 7 kHz strength disorder (c). d) shows a cut along z (integrated along x) with 
a bi- modal distribution fit by a single Gaussian (i), with a Bose function fit to the 
thermal portion (ii), and with a Gaussian fit to the central peak with the thermal 
portion removed (iii). A well-separated vortex originally present in one of the planes 
is indicated by the dislocation of the phase of the fringe pattern near the center of the 
cloud in (a). 

would like to get information about the properties (particularly phase coherence) of the 
trapped atoms by looking at the atoms after 22 ms of expansion. 

To disentangle the effects of density and phase modifications in the trapped gas 
and their effects on our measurements, we have first verified that the disorder has no 
effect on the temperature (see Supplementary Material). If only the phase coherence 
of the trapped samples becomes disordered, we expect that the integrated visibility 
of the interference pattern observed after expansion will decrease. A smaller local 
trapped density (induced by a local minimum in the disorder) will also decrease the 
local visibility after expansion, but if those missing atoms remain in the central peak 
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Figure 3. Disorder-Phase Space Phase Diagrams Color maps showing integrated 
visibility over the central portion of the interference pattern (a) and coherent fraction 
(b) in our images as a function of disorder strength and phase space density. The color 
maps are qualitatively similar but the drop in coherent fraction at high disorder and 
high D (top right corner of the color maps) is much less than the corresponding drop 
in visibility. 



during the expansion, the integrated visibility of the central peak will not change. The 
disorder could also disrupt overall coherence such that atoms which would have appeared 
in the central peak after expansion are redistributed to the wings of the distribution 
during the expansion. We can measure this effect by counting the fraction of atoms in 
the central peak in our images. 

3. Analysis 

Each image of the expanded cloud is analyzed to obtain coherent fraction, temperature, 
total atom number, and integrated fringe visibility. Total atom number is obtained by 
integrating the density distribution in each image. To obtain temperature and coherent 
atom number, we use the observed two-component distribution along the z-axis in the 
images. We integrate along the x-direction and then exclude the central portion of 
the cloud, fitting the remaining distribution along z to a Bose function. The width 
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of this Bose function is used to obtain the temperature. We subtract off the Bose 
fit, assuming that the thermal component in the x-direction has the same width as 
the coherent component, and then integrate the remaining atom density to obtain the 
coherent atom number. The integrated fringe visibility is determined by taking the ID 
Fourier transform of the coherent atoms along x and finding the amplitude A(z) and 
phase (p(z) of the Fourier peak at the frequency of the fringes. 

Figure |3k) shows the average fringe visibility J A{z)e l ^ z 'dz integrated over 120% 
of the Gaussian width of the central portion of each cloud at different peak phase 
space densities and disorder strengths. The D shown here is what the peak D 
in each of our planes would be at the measured temperature and atom number in 
the absence of disorder, calculated with a local density approximation (LDA, see 
supplementary material). Since the local properties of the (non-disordered) gas are 
completely determined by D, this is the appropriate quantity to use along one axis. We 
have also measured the fraction of atoms in the coherent central part of the measured 
bi-modal distribution (Fig. [3b). We observe a similar phase diagram for this central 
coherent fraction as for the integrated visibility. However, the effect of the disorder is 
not as strong, most notably in its effect at high D and high Af/h. 

We can quantify the difference in fractional decrease in visibility vs. fractional 
decrease in coherent fraction as a function of disorder strength (Fig. H]). Ideally, we 
would like to compare the asymptotic values that the visibility and coherent fraction 
approach at high D and plot them as a function of disorder. However, it is unclear if 
either visibility or coherent fraction is approaching a constant value at high D, even 
in the absence of disorder. Instead, we average the fractional change across all phase 
space densities at each disorder strength. This serves to weaken the signal, as we expect 
the role of disorder to decrease at low D. Still, we find that once the disorder strength 
exceeds 4 kHz, the fractional decrease in visibility is ~30% more than the fractional 
decrease in coherent fraction. 

In regions of the sample where D is high enough to be at least quasi-condensed, we 
can describe the system with an order parameter having both an amplitude and a phase 
- these are the atoms in the central peak of the expanded distribution. An explanation of 
Fig. H] is that the disorder has affected the phase of this order parameter in each system 
while not affecting the amplitude of the order parameter nearly as much, as predicted 
by the dirty Boson model [H [5j E]. This theory shows that disorder generates phase 
fluctuations in a low temperature superfluid that drive the system into an insulating 
state. In a quasi-2D ultracold gas, it is expected (but not yet experimentally verified) 
that the gas is locally superfluid when D > D c , which may make this model applicable. 
Until now, the effect of disorder on quasi-2D superfluids described by this model has 
been limited to liquid 4 He thin films [29] or thin film superconductors [30], with no 
access to the phase of the superfluid order parameter. In our system, we can measure 
the phase of the order parameter, but we do not observe a clear indication of a transition 
at high disorder strengths. Instead, we see the visibility level off to a D-dependent value 
at high disorder strength. 
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Figure 4. Fractional Change Due to Disorder Average fractional change for the 
coherent fraction and integrated visibility as a function of disorder strength. Each 
data point is the average over D at that disorder strength, and error bars are 90% 
confidence intervals for that average. 



This leveling may be related to confinement or finite-sized system effects. The 
relevant length scales are the system size (20 x 90 /jjo), the disorder length scale (~0.53 
/mi), and the healing length (0.25 /mi). We can investigate the effect of larger disorder 
correlation lengths by decreasing the numerical aperture of our disorder optical system. 
Models based on random Josephson junctions (31] or percolation [32j may also help 
explain these results, but a clear theoretical understanding is made difficult by the 
non-uniform system density and the exact nature of the coherence in the central peak. 
Other ultracold atom disorder experiments may have also observed similar decoupling 
of phase and condensate fraction [331 EH EH], but there are significant differences 
in the dimensionality and underlying structure in these experiments, making direct 
comparisons difficult. 

We observe well-separated vortices in some of our images. The statistics of these 
vortices may be used to further characterize the effect of disorder in 2D systems, 
possibly indicating pinned vortices, as in type II superconductors [35]. Preliminary 
measurements suggest that the percentage of realizations with one or more vortices 
present does not depend on the amount of disorder, but more data must be taken to 
support any conclusions. We can also translate the disorder relative to the confining 
potential, allowing for the exploration of vortex pinning and transport. We have shown 
that disorder can drive phase fluctuations in quasi-2D systems, a direct observation 
of the mechanism predicted for the disorder- driven SIT in the dirty Boson model. 
This has relevance to systems such as superconducting thin films and high-T c cuprate 
superconductors, which may undergo phase transitions by this mechanism. This result 
may bring us closer to understanding the mechanism behind phase transitions in high-T c 
sup er conduct ors . 
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Appendix 

Appendix A.l. Disorder Calibration. 

The correlation length of the disorder was calibrated by imaging the disordered intensity 
pattern J(r) using a microscope objective and CCD camera. We computed the 
autocorrelation function C(r) = (/(r )/(r + r)), where the brackets represent an 
average over ro, and used the Gaussina width of the function to determine the correlation 
lengths. For the transverse autocorrelation length, we used a single image of the beam 
at the focus of the final lens. To obtain the axial autocorrelation length, the imaging 
system was mounted on a translation stage with ~30 nm step size, and a series of images 
was taken and then corrected for transverse jitter before the autocorrelation function of 
a slice along the axis was computed. 

To calibrate the strength of the disorder, we pulsed the disordered potential 
on a thermal cloud of atoms and observed the increase in momentum distribution 
width. If the pulse is short compared to atomic motion timescales, it gives an impulse 
AP(r) = V£/(r), where U(r) is the potential strength of the disorder, proportional to 
the intensity of the beam. We modeled the process using a Monte Carlo simulation 
and compared the actual increase in momentum distribution width to our simulation in 
order to obtain disorder strength as a function of power in the disordering beam. 

Appendix A. 2. Color Maps 

The bin sizes for the color maps were chosen so that almost all bins had 2 or more points, 
with some bins in the middle of the maps containing as many as 50 points. A total of 
1300 images were analyzed to produce the color maps. Each bin is colored according to 
a Gaussian weighted average over the entire data set, with the Gaussian width set to 
1/3 of the bin size. 

Appendix A. 3. Disorder Effect on Temperature Measurement 

We took several hundred alternating images with our lattice and disorder to determine 
the effect of the lattice and disorder on the measured temperature of the thermal cloud. 
The first image was taken with our typical procedure described above, while the second 
image in each set followed the same timing, but the lattice and disordering beams were 
never turned on. For the 3D BEC image, the center was masked and we fit ID cuts 
of the thermal background in two directions to Gaussians to extract the temperature. 
We show in Fig. IA1I the ratio of the temperature measured in concurrent shots for a 
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Figure Al. Effect of Lattice and Disorder on Temperature Ratio of 
temperature measurements in concurrent shots with and without the lattice and 
disorder. Each point is the average of ~ 5 shots at the same value, with error bars 
the standard deviation. For a given temperature set, the quoted temperature is the 
average temperature of the shots without the lattice and disorder, with the standard 
deviation giving the error. The flatness of the curves indicates that increasing disorder 
has no effect on the measured temperature of the sample. 



variety of disorder strengths. The ratio looks flat as a function of disorder, indicating 
that higher disorder strengths have a negligible effect on our temperature measurement. 
We repeated the measurement at various average temperatures to be sure that there 
was no temperature dependence to the curve. A combination of adiabatic heating and 
spontaneous emission heating increase the temperature in the lattice and disorder by a 
factor of ~ 1.4 from its 3D BEC value. 

Appendix A. 4- Phase Space Density and Chemical Potential 

Once each of our images is analyzed to obtain temperature and total atom number, these 
quantities (along with the properties of the trapping potential) can be used to calculate 
phase space density and chemical potential in the analogous disorder-free system. It 
is easiest to divide quasi-2D systems into three distinct phases: normal, fluctuation, 
and Thomas-Fermi regimes, corresponding to low, medium, and high densities. In the 
normal regime, the density can be obtained using the following procedure from Refs. 
[201123]. 

The 2D in-plane density n(r) in each vibrational level j (assumed to be quantum 
harmonic oscillator eigenstates) along x is 

%(r0 = -In (l-e^ k * T )/\ 2 dB , (A.l) 

with 

1 fAixTi 2 \ 

^(f) = ^ai _ _ m(u y + ^2) _ jhUx _ £ 2 i^— -a s f m (r) J, (A.2) 
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where kg is Boltzmann's constant, Xjb = y2nh /(raksT) is the thermal deBroglie 
wavelength, ^p lohal is the peak chemical potential, m is the mass of 87 Rb, w» is the 
harmonic trapping frequency in the given direction, h is Planck's constant, and a s is 
the s-wave scattering length of 87 Rb. In Eq. IA.21 fji is the normalized density overlap 
integral between harmonic oscillator eigenstates j and /. Eq. IA.1I can be solved via an 
iterated fixed point method, and converges after a few iterations. Note that phase space 
density Dj(r) = nj{r)X 2 dB . 

In the high-density Thomas-Fermi region of the cloud, the density is given by Eq. 
IA.ll except that the factor of 2 in front of the sum in Eq. IA.2I is not there. That factor 
of 2 is there in the normal phase because density fluctuations are important, so that 
(n 2 ) = 2 (n) 2 . When the density gets very high, density fluctuations are suppressed, 
and (n 2 ) = (n) 2 , removing the factor of 2 in the sum term of Eq. IA.2I 

In between these two densities, the relationship is more complicated. In the strictly 
2D case, where only the lowest harmonic oscillator eigenstate is occupied, the phase 
space density in this fluctuation regime is calculated with a Monte Carlo method, and 
the results are tabulated by Prokof'ev and Svistunov [36]. These results were extended 
to the quasi-2D case for finite occupation of excited vibrational modes by Holzmann, 
Chevallier, and Krauth [37]. They found that one could define a correlation density AD 
- the difference between the true phase space density in the fluctuation regime and what 
is predicted by Eq. IA.1I The correlation density is a function only of the mean field gap 

A mf = -In (l - e- Do ) , (A.3) 

with D Q the ground vibrational state phase density calculated using Eq. IA.ll and is 
given with 10% accuracy up to the BKT transition density D c to be 

AD ~ - ( -1 + J—\ ] 9 , - . (A.4) 

We extended this approximation D(r) = D m f(r) + AD(f) [D m f from Eq. lA.lj) through 

the fluctuation region to densities higher than D c , until it smoothly matched up to the 

Thomas-Fermi approximation. 

The method detailed above allows us to self-consistently calculate the density over 

any trapped sample. Since we measure total atom number in both of our planes, we find 

a value of ^ lobal which gives us (at the measured temperature) a calculated integrated 

density matching half of the total measured atom number. The peak phase space density 

shown in the figures in the main text is the calculated value of D(0). 
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